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Abstract. We present a new methodology to determine the expansion history 
of the Universe analyzing the spectral properties of early-type galaxies (ETG), 
based on the study of the redshift dependence of the 4000 A break. In this 
paper we describe the method, explore its robustness using theoretical synthetic 
stellar population models, and apply it using a SDSS sample of ~14 000 
ETGs. Our motivation to look for a new technique has been to minimize 
the dependence of the cosmic chronometer method on systematic errors. In 
particular, as a test of our method, we derive the value of the Hubble constant 
H = 72.6 ± 2.9(stat) ± 2.3{syst) kmMpc"^- 1 (68% confidence), which is not 
only fully compatible with the value derived from the Hubble key project, 
but also with a comparable error budget. Using the SDSS, we also derive, 
assuming w = constant, a value for the dark energy equation of state parameter 
w = -1 ± 0.2(stat) ± 0.3(syst). Given the fact that the SDSS ETG sample only 
reaches z ~ 0.3, this result shows the potential of the method. In future papers 
we will present results using the high-redshift universe, to yield a determination 
of H (z) up to z ~ 1. 



1. Introduction 

In the last decades many efforts have been done in trying to unveil and understand the 
expansion history of our Universe. Many different and complementary observations, 
such as direct supernova measurements of the deceleration parameter (Riess et al. 
1998, Perlmutter et al. 1999) as well as indirect measurements based upon a 
combination of results from the cosmic microwave background (CMB) (Jungman et al. 
1996, Miller et al. 1999, De Bernardis et al. 2000, Hanany et al. 2000, Jaffe et al. 2001, 
Halverson et al. 2002, Mason et al. 2003, Benoit et al. 2003, Goldstein et al. 2003, 
Spergel et al. 2003, Spergel et al. 2007, Reichardt et al. 2009, Dunkley et al. 2009), 
large-scale structure (LSS) (Percival et al. 2001, Dodelson et al. 2002), and the Hubble 
constant (Freedman et al. 2001), pointed out that at present time the expansion of the 
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Universe is accelerating. This could be explained either by considering that gravity at 
some large scales cannot be described by the standard general relativity and/or that 
the Universe is filled with some sort of negative-pressure "dark energy" that drives 
the accelerated expansion (Peebles et al. 2003, Padmanabhan 2003, Copeland et al. 
2006, Fricman et al. 2008, Linder 2008, Caldwell et al. 2009); either way, it requires 
new physics beyond general relativity and the standard model of particle physics. 

The simplest possibility is to extend Einstein's equation with a cosmological 
constant, or equivalently, to hypothesize a fluid with an equation-of-state parameter 
w = p/p = — 1 (with p and p the pressure and energy density, respectively). However, 
it may well be that the cosmological "constant" actually evolves with time, in which 
case ic/-l, and there are a variety of possibilities to believe that this might be the 
case (Ratra et al. 1988, Caldwell et al. 1998). Precise measurement of w(z) (with, 
in general, a parameterized redshift dependence) or, relatedly, the cosmic expansion 
history H(z), has thus become a central goal of physical cosmology (Peacock et al. 

2006, Albrecht et al. 2006). 

At present, the most promising techniques to determine the cosmic expansion 
history are supernova searches (Riess et al. 1998, Perlmutter et al. 1999), baryon 
acoustic oscillations (BAO) (Seo et al. 2003, Eisenstein et al. 2005, Percival et al. 

2007, Pritchard et al. 2007), weak lensing (Refregier et al. 2003), and galaxy clusters 
(Haiman et al. 2001). All these techniques have different strengths and suffer from 
different systematics and weaknesses. As argued in the ESO/ESA and Dark Energy 
Task Force reports (Peacock et al. 2006, Albrecht et al. 2006), robust conclusions 
about the cosmic expansion history will likely require several independent methods to 
allow for cross checks and control of systematics. 

As suggested by Jimenez & Loeb (2002), an interesting and complementary 
approach to this problem is the study of the change in the age of the Universe as 
a function of redshift. The strength of this method is that it avoids the common 
weakness of the other techniques, i.e. the reliance on integrated quantities to determine 
the expansion history. However, to fully exploit this approach, it is required to find 
a population of standard clocks able to trace the evolution of the relative age of the 
Universe. 

Many recent works gave growing observational evidence that the most massive galaxies 
contain the oldest stellar populations up to redshifts of z ~ 1 — 2 (Dunlop et al. 1996, 
Spinrad et al. 1997, Cowie et al. 1999, Heavens et al. 2004, Thomas et al. 2005, 
Cimatti et al. 2008, Thomas et al. 2010) and have less than 1% of their present 
stellar mass formed at z < 1. Studying the mass function of early- type galaxies, it 
has been found that its massive end evolves very mildly (or not at all) from z ~ 0.7 
to z — (Cimatti et al. 2006, Pozzetti et al. 2010). This population of massive, red, 
passive early-type galaxies (ETGs) consists of the oldest objects in the Universe at 
each redshift. The differential ages of these galaxies should thus be a good indicator 
for the rate of change of the age of the Universe as a function of redshift up to z ~ 2, 
given that it is commonly accepted that the bulk of their stars was formed at z > 2 — 3. 
There are some works that try to constrain the cosmic expansion history using this 
approach (Ferreras et al. 2001, Jimenez et al. 2003, Krauss et al. 2003, Capozziello 
et al. 2004, Stern et al. 2010, Crawford et al. 2010). However, a problem of this 
technique is that it relies on the determination of a non-observable parameter using 
spectra and/or spectral energy distributions, i.e. the age of the stellar population, 
which presents strong and well known degeneracies with other parameters that are 
evaluated simultaneously, such as the metallicity, the star formation history, the dust 
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content. 

The aim of this project is to search for features in the spectra of ETGs that 
minimize the impact of systematics on the estimate of the relative ages of galaxies, 
which is the relevant quantity to determine H(z) using the cosmic chronometer 
method. As we will demonstrate below, the 4000 A break is a feature correlated almost 
linearly with age at fixed metallicity. The 4000 A break (D4000) is a spectral feature 
produced by the blending of a large number of absorption lines in a narrow wavelength 
region, that creates a break in the continuum spectrum at 4000 A restframe. The main 
contribution to the opacity comes from ionized metals. In hot stars, the metals are 
multiply ionized and the opacity decreases, so the D4000 will be small for young 
stellar population and large for old, metal rich galaxies. A measure break amplitude 
was firstly introduced as the ratio of the average of the flux densities below and above 
4000 A (Hamilton, 1985); he was also the first to suggest that this feature can be used 
as an age indicator for a galaxy. Since then, different index definition were proposed, 
depending on the width of the red and blue regions where the fluxes are evaluated, 
and many studies have been done on this feature. Balogh et al. (1999) introduced a 
definition of D4000 with narrower bands, supporting this choice with a lower sensitivity 
to reddening effects. Bruzual & Chariot (2003) created a model for computing the 
spectral evolution of stellar populations (BC03), capable also to trace the dependence 
of D4000 on age and metallicity as a function of age, reddening and star formation 
history. Kauffmann et al. (2003) found, from SDSS observations, that the distribution 
of D4000 is strongly bimodal, showing a clear division between galaxies dominated by 
old stellar population and galaxies with more recent star formation. 

In this paper we analyze almost 14 000 ETGs extracted from the SDSS survey, 
in the redshift range 0.15 < z < 0.3. We will show that for these old, red and passive 
objects the D4000 is, if the metallicity is known, an optimal indicator that correlates 
almost linearly with age, depending only marginally on the adopted synthetic stellar 
population models or star formation history (SFH). It is therefore possible to use 
the differential D4000 evolution of the galaxies as a tracer for their differential age 
evolution, and thus to set constraints to cosmological parameters in an almost model- 
independent way. The paper is organized in the following way. In Sect. [5J we present 
the data sample, describe the different criteria with which ETGs have been selected, 
and give information about the D4000 definition adopted; we also discuss the mass 
and metallicity evaluation, and give more details about the distribution of various 
parameters that enhance the robustness of our selection. In Sect. G2 we describe the 
analysis done on the D4000-age relation obtained both using Bruzual & Chariot (2003) 
and Maraston & Stromback (2010) stellar population synthesis models, showing that 
in the range of D4000 probed by our data the linear approximation works extremely 
well and evaluating the slope of this relation as a function of metallicity. We therefore 
obtain the theoretical relation linking the D4000 with the redshift as a function of the 
various cosmological parameters. In Sect. 2] we show the D4000-Z relations obtained 
for different mass subsamples of our data, obtaining a clear redshift evolution for all the 
sample and a strong evidence of mass-downsizing. In Sect. [5]we show the results of the 
evaluation of Ho and of a constant w, comparing them with the recent determination 
done by Riess et al. (2009), by the WMAP 7- years analysis and by Amanullah et al. 
(2010). We also discuss the robustness of our results against possible systematics from 
which our analysis may suffer, such as the dependence on the estimate of the absolute 
age of a galaxy from the D4000„, on the star formation histories considered, on the 
stellar population synthesis models assumed, on the metallicity evaluation and on the 



Cosmological constraints with spectroscopic evolution of cosmic chronometers 4 

not perfect assumption of the passive evolution of our ETGs, reporting at the end our 
estimates of H and w both quoting statistical and systematic errors. 

2. The sample 

2.1. Early-type galaxies selection 

The sample of galaxies analyzed in this paper has been selected from the spectroscopic 
catalog of the Sixth Data Release of the Sloan Digital Sky Survey (SDSS-DR6). This 
survey provides u, g, r, i, z photometry and spectra for almost 800 000 galaxies over 
~ 7000 square degrees. Galaxies have been extracted to have a Petrosian magnitude 
r < 17.77. The spectra are taken using 3-arcsec diameter fibres, positioned as close 
as possible to the centers of the target galaxies. The flux (and wavelenght) calibrated 
spectra cover the range 3800-9200 A, with a resolution of R ~ 1800. To obtain a wider 
wavelength coverage for the photometry, we use a match between SDSS-DR6 galaxies 
and 2MASS survey provided to us by Jarle Brinchmann (private communication), 
with which we extend the SDSS photometry with the J, H and K bands. In this way 
we obtain for each galaxy a 8-bands spectral energy distribution (SED) from u to K, 
very important for a more robust determination of mass and of colors of the galaxy. 
As explained before, it is fundamental to have a population of standard clocks in 
order to do precision cosmology with the differential age technique. In a previous 
paper (Moresco et al. 2010), we checked on the zCOSMOS 10k bright sample that a 
simple color selection is not restrictive enough to exclude starforming galaxies from 
the early- type galaxies selection. Following the analysis explained in that paper, we 
therefore decided to apply a more restrictive selection, on the basis of photometric 
and spectroscopic informations, to obtain a uniform, passive population of early-type 
galaxies (ETGs); we selected from SDSS-DR6 all galaxies matching these criteria: 

• no strong emission lines: we decided to select only galaxies with 
£W ( [O/I] A3727) < 5 A and EW (H a ) < 5 A. These two spectroscopic fea- 
tures are strongly correlated with ongoing or recent episodes of star formation 
history. This upper limit in the equivalent width allows us to select galaxies that 
have experienced a limited (or null) amount of star formation in the last 1-2 Gyrs, 
thus choosing the most passive sample. 

• early-type SED: following the procedure explained in Zucca et al. (2006), we 
performed a best-fit modeling of the multi-band photometry of these galaxies 
using the empirical set of 62 SEDs described in Ilbert et al. (2006), where these 
SEDs were derived by interpolating between the four local observed spectra of 
Coleman et al. (1980) (from the old stellar population in M31 and M81 to Sbc, 
Scd, and Im SEDs) and two starburst SEDs from Kinney et al. (1996). In this 
way we selected only those galaxies that match a early-type SED template. 

For each galaxy, the 4000 A break has been taken from the MPA-JHU DR7 release 
of spectrum measurements \http:/ /www. mpa-g arching. mpg. de/SDS S/DR7/) , which 
provides a complete emission line analysis for the SDSS Data Release 7 (DR7). We 
decided to use the D4000 definition proposed by Balogh et al. (1999), that uses 
narrower bands in order to be less sensitive to reddening effects (3850-3950 A and 
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Figure 1. Mean stacked spectrum of 10 ETGs randomly selected from the 
extreme tail of the D4000 n distribution. The red shaded areas show the ranges 
used to evaluate the D4000 n and the red horizontal lines show the mean of the 
flux. 



4000-4100 A); we denote this index as D4000 n : 

F red (\% ue -\?™)jSF u dA 

m000n = -p± = — - — -— — (1) 

tblue (\red _ ^ed) J^y F^dX 

In Fig. Q]we plot the average stacked spectrum of 10 random galaxies, selected to 
be in the highest tail of the D4000„ distribution (D4000„ > 2.1) and for which the 
-D4000„ has been measured with high accuracy (S/N ~ 70). The dispersion of the 
-D4000„ of the galaxies used for creating the stacked spectrum is < 0.01. In the figure 
are shown also the ranges used to evaluate the -D4000 n and the mean fluxes in the 
corresponding regions. 

As pointed out by Bernardi et al. (2006), the spectrophotometric calibration of the 
SDSS spectra is less reliable around 4000 A, due to the throughput of the spectrometer. 
This means that the value of the D4000„ can be affected by this problem for galaxies 
in the lowest redshift range (Bernardi, private communication). In order to avoid this 
potential source of uncertainty, we decided to select our final sample at z > 0.15. 
Moreover, since the number of galaxies decreases very rapidly at z > 0.3, we excluded 
galaxies beyond z=0.3. Thus, the final redshift range of our sample is 0.15 < z < 0.3. 
As explained in the forthcoming section, the metallicity estimate for these galaxies has 
been taken from the SDSS-DR4 analysis performed by Gallazzi et al. (2005). After 
the match between our sample and the galaxies for which the metallicity estimate is 
publicly available, we end up with a sample of 13 987 galaxies, for each one having 
a wide photometry coverage (u to K), spectroscopic redshift, spectral index measures 
(£>4000„, EWo(Ha) and EW ([OII]X3727)), metallicity and mass estimates. 



Cosmological constraints with spectroscopic evolution of cosmic chronometers 



6 




Figure 2. In the upper panels, distributions of r values (left panel) and 
mctallicity-rcdshift relations (right panel) in different mass subsamples. In the 
lower panels, distributions of the EWo(Ha) (left panel) and EW ([OII] A3727) 
(right panel). The red lines represent the median of the distributions. 



2.2. Mass, metallicity and star formation history estimation 

To study the properties of our sample of ETGs, we used both photometric and 
spectroscopic information. To retrieve information about the metallicities of those 
galaxies, we used the analysis performed by Gallazzi et al. (2005), in which 
they derived stellar metallicities for a sample of SDSS-DR2, afterward extended 
to SDSS-DR4 (www.mpa-garching.mpg.de/SDSS/DR4/Data/stellarmet.html). Their 
constraints are set by the simultaneous fit of five spectral absorption features (D4000, 
Hj3 and HS a + H^ a as age-sensitive indices and [Mg2Fe] and [MgFe]' as metal- 
sensitive indices, all of which depend negligibly on the a/Fe ratio), which are well 
reproduced by Chariot & Bruzual (2003, BC03) population synthesis models (see 
Gallazzi et al. 2005 for further details). 
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Figure 3. Distributions of the difference (ZQ a u azzi — Zvespa) (left panel) and 
of the percentage difference (Z G allazzi - Z v E s pa) /ZGallazzi- 



To verify that the metallicity estimation is not biased by the choice of the particular 
star formation history done by Gallazzi et al. (2005), we decided to compare their 
estimates with the ones done with VESPA code (Tojeiro et al. 2007). The approach 
of VESPA is to recover robust star formation and metallicities histories using the 
full spectral range of a galaxy from synthetic models; in order to estimate those 
quantities, it allows for a completely free-form star formation, imposing no prior on 
it. This code has been applied to the entire SDSS-DR7, obtaining a catalog of stellar 
masses, detailed star formation and metallicity histories (Tojeiro et al. 2009, www- 
wfau. roe. ac. uk/vespaf) . 

Despite the completely different assumption on the star formation histories, 
the distribution of the difference (ZQ a u azzi — Zyespa) is a gaussian with mean 
value -0.01 and dispersion 0.1; the distribution of the percentage difference 
(ZGallazzi - Z v E s pa) / Zcaiiazzi has a mean value of 4.8%, with a 51.3% dispersion. So 
there exists no bias between the two metallicity evaluations, and this analysis ensures 
the reliability of Gallazzi et al. (2005) metallicity estimates (see Fig. [3]). 

The masses are estimated by performing a best-fit to the multicolor spectral 
energy distribution (SED), using the observed magnitudes in 8 photometric bands 
from u to K; to fit our data, we decided to use a grid of models from BC03 stellar 
population models. The grid has been created with delayed exponential SFH with 
SFR(t,r) oc (t / T 2 )exp(—t / t) (where r is chosen in the range 0.05 < r < 1 Gyrs), 
reddening in the range < Ay < 1 and free ages (see Pozzetti et al. 2010 for 
more details). Many recent works on early-type galaxies evolution has shown clearly 
the existence of a mass-downsizing effect, i.e. more massive galaxies have assembled 
their mass before less massive ones. Cimatti et al. (2006) analyzed a sample of 
COMBO-17 and DEEP2 ETGs, finding that the amount of evolution for the ETG 
population depends critically on the range of luminosity and masses considered. 
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Scarlata et al. (2007b) show similar results for the photometric survey COSMOS, using 
both morphologically and photometrically selected subsamples of early-type galaxies. 
Pozzetti et al. (2010), studying galaxies obtained by the zCOSMOS survey confirm a 
mass dependent evolution of galaxies up to z ~ 1 in particular for ETGs, i.e. galaxy 
number density increases with cosmic time faster for less massive galaxies and < 20% 
for massive ones. Ilbert et al. (2010), extend the result at z > 1, selecting ETGs in 
SCOSMOS (Sanders et al. 2007) up to z p hoto = 2, finding a more rapid evolution 
for massive galaxies at higher redshifts (by a factor of 15-20 between z = 1.5 — 2 and 
z = 0.8 — 1). The analysis performed by Thomas et al. (2005, 2010) on a sample taken 
from the SDSS survey reveals that galaxies above log(M /Mq) ~ 12 have completed to 
assemble their mass already at z w 3, while this process happens only around z 1 for 
galaxies with log(M /Mq) ~ 10. Moresco et al. (2010) confirmed this result analyzing 
a sample of ETGs extracted from the zCOSMOS 10k bright spectroscopic sample (see 
Lilly et al. 2009 for further details on the parent sample); they found that the mean 
redshift of formation for galaxies with log(M /Mq) ~ 11 is around z » 2, while for 
galaxies with log(M /Mq) ~ 10 it is around 

In order to select the oldest population from our ETGs sample, we therefore decided 
to select the high-mass tail of the mass distribution, considering log(M /Mq) > 10.75. 
Moreover, to avoid possible bias from mass-downsizing effect, we decided to split 
our sample in narrow mass bins, to have a more homogenous sampling of redshift of 
formation and metallicity for these galaxies. We divided our sample in four mass bins, 
with 10.75 < log{M/M & ) < 11, 11 < log(M/M Q ) < 11.25, 11.25 < log(M/M Q ) < 
11.5, and log{M/Mo) > 11.5. The division our ETGs into four subsamples of the 
width of Alog(M /Mq) = 0.25 has been done to keep the statistical significance of 
the subsamples high enough, having for each subsample N ga i > 3500 (except for the 
highest mass bin, for which N ga i ~ 1500); we adopted the cut log{M / Mq) > 10.75 in 
order both to select the most massive galaxies of our sample and to have four mass 
bins with a sufficient number of galaxies to perform a reasonable statistical analysis. 
We checked that the difference of the median mass along the redshift range is 0.11 
dex on average and that the median metallicity is constant in the considered range 
redshift within a 5% percent level, on average (see upper-right plot of Fig. 

The t parameter of the SFH obtained from the best-fit to the SEDs of the galaxies 
shows a distribution characteristic of passive galaxies, with median values below 0.2 
Gyrs for all the mass subsamples (see upper-left plot of Fig. [2]). This result is in 
agreement with many ETGs analysis (Cimatti et al. 2008, Longhctti et al. 2008, 
Gobat et al. 2008), for which the majority of massive field early- type galaxies should 
have formed their stellar content around z > 2 over short (i.e. r < 0.1 — 0.3 Gyrs) star 



mass range 


median mass 


median Z/Zq 


median EWo(Ha) 


median EWo{[OII]) 


# galaxies 


log(M/M Q ) 


log(M/M e ) 


Z/Zq 


[A] 


[A] 




10.75-11 


10.9 ±0.02 


1.061 ±0.006 


0.463 ± 0.009 


1.24 ±0.03 


3452 


11-11.25 


11.12 ±0.03 


1.088 ± 0.004 


0.376 ± 0.007 


1 ± 0.02 


5429 


11.25-11.5 


11.25 ±0.02 


1.1 ±0.006 


0.305 ±0.007 


0.81 ±0.02 


3591 


> 11.5 


11.63 ±0.04 


1.144 ±0.009 


0.27 ±0.01 


0.78 ±0.03 


1515 



Table 1. Median value of the mass (in logarithmic units), of the metallicity (in 
solar metallicity units), and of the Ha and [OI I]A3727 equivalent widths of the 
ETGs in different mass subsamples. We also quote for each bin the number of 
galaxies. 
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formation time-scales. We also studied the distributions of the Star Formation Rate 
(SFR) obtained from the SED-fitting, and found a median SFR equal to zero M Q /yr 
for all the mass subsamples, with a scatter of 0.11, 0.17, 0.18, 0.2 M Q /yr, respectively 
for the four mass subsamples. 

Since our analysis depends strongly on the selection of a sample of galaxies 
passively evolving, we analyzed in detail also the emission lines distributions, because 
they can be linked with star formation or AGN activity. We find that the median 
values of those distributions are EW (Ha) < 0.5 A and EW a ([OII}\3727) < 1 A, 
with tails up to 2 A for the EW (Ha) and up to 4 A for the EW {[0 1 I}\37 '27) (see 
Tab. [IJ. The median value of the specific star formation rate, obtained from the 
SED-fitting, for the galaxies with a significant detection of EWo(Ha) is however null. 
Moreover, we verified that the Z?4000„ is not significantly different between galaxies 
with EWo(Ha) < and EW (Ha) > 0, finding respectively L>4000„ = 1.948 ± 0.002 
and L>4000„ = 1.944 ± 0.001. 

We also checked the distribution of the eClass spectral parameter, that is a SDSS 
parameter quantifying the activity of a galaxy using a PCA analysis of the spectra of 
the galaxies; it goes from -0.35 to 0.5 for early- to late-type galaxies. For each mass 
subsample, it shows a distribution always below -0.1, so our selection results even 
stricter than the one done from Bernardi et al. (2006), in which they chose ETGs 
with eClass < 0. In Fig. [5] we show in the upper panels the r distribution (left 
panel) and the metallicity-redshift relations (right panel), while in the lower panels 
the distribution of the EW (Ha) (left panel) and EW Q ([OII}\3727) (right panel), 
each one for the four samples; the red lines represent the medians of the distributions. 
In Table [T] are reported the median values of mass, metallicity and equivalent widths 
of the emission lines for the different mass subsamples, with their errors; we report 
also in the same table the number of galaxies in each bin. 

Another possible systematic trend with redshift that could, in principle, bias our 
analysis is the effect of using in SDSS a fixed 3 arcsec fibre, which maps different 
physical sizes depending on the observed redshift. However Gallazzi et al. (2005) 
analyzed in details the effect of this aperture bias, finding no significant difference in 
metallicity between high- and low-concentration galaxies (see their Sect. 3.2). We also 
checked that the difference in D4000 between high- and low-concentration galaxies is 
always < 2%. 

3. Linking the 4000 A break to the expansion history of the Universe 

3.1. Calibration of the D4000 n -age relation 

In general the D4000„ is an index that is strongly sensitive both to metallicity, star 
formation history and age of a stellar population. In order to break this degeneracy, 
we created a library of Z?4000„ as a function of age using stellar population synthesis 
models with different metallicity, star formation history (delayed exponential SFR 
with r = 0.05, 0.1, 0.2, 0.3 Gyrs), and ages. In order to study the possible systematics 
due to the choice of a particular model, we created two different libraries, one using 
Bruzual and Chariot (2003) stellar population synthesis models (hereafter BC03) and 
one using the new Maraston & Stromback (2010) stellar population synthesis models 
(hereafter Mastro), which are based on the MILES models (Sanchez-Blazquez et al. 
2006). The metallicities studied with the two models are Z/Z Q = 0.4, 1, 2.5 for BC03 
and Z/Zq = 0.5,1,2 for Mastro. The resolution of the two models is similar, 3 A 
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Figure 4. D4000 n -age relation for BC03 high resolution models (left panel) 
and M astro models (right panel). In the upper panels the colored lines represent 
models with different metallicities, with oversolar metallicity in red (respectively 
Z/Zq = 2.5 and Z/Zq = 2), solar in green, and undersolar in blue (respectively 
Z/Zq = 0.4 and Z/Zq = 0.5). The blue and the red shaded area represent 
the ranges of D4000n in the models (respectively the low D4000 n and the high 
Z)4000 n regime). The lower panel show a zoom of the interested area, where 
the models are shown in gray and the colored lines are the fit to the models. 
The dotted line shows where the change between the low D4000 n and the high 
D4000n regime. It is possible to see the minor dependence of the slopes of the 
D4000 n -age relation on the SFH assumed both by visually comparing the average 
slopes (colored lined of the lower panels) and the slope of each single model (grey 
lines of the lower panels) and by looking the small errors associated to the mean 
values (see Table [2{. 

across the wavelength range from 3200 A to 9500 A for BC03, and 2.3 A across 
the wavelength range 3525 A to 7500 A for Mastro; moreover the resolution of both 
models is comparable to the resolution reachable with the SDSS spectrograph, which 
has R ~ 1800 between 3900-9100 A. 

The choice of the grid of star formation histories and metallicities is motivated 
by the fact that we want to analyze ETGs from the SDSS survey: a lower value in 
the range of metallicity would probe ZM000„ values that are much lower than the 
average values found in our galaxy sample, and the SED-fitting analysis of our ETGs 
shows that the majority of the sample is best-fitted with SFHs with low values of r. 
Moreover, the choice of a passive sample strongly supports our decision of using these 
low t values. 

The upper panel of Fig. [4] shows the BC03 and Mastro £>4000„-age relations; 
the models with the higher metallicity (Z/Zq — 2.5 and Z/Zq = 2 respectively) are 
in red, the ones with Z/Zq = 1 in green, and the ones with the lower metallicity 
(Z/Zq — 0.4 and Z/Zq = 0.5 respectively) in red. For each metallicity, the models 
are plotted with a continuous line for r = 0.05 Gyrs, with a dotted line for r = 0.1 
Gyrs, with a dashed line for r = 0.2 Gyrs and with a long-dashed line for r = 0.3 
Gyrs. 
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We decided to study these relations in the D4000„ range spanned by our data, 
roughly 1.8 < D4000„ < 2. For both the models considered, we found that 
the Z?4000 n — age relation presents, at each metallicity, two different slopes: one 
characteristic of the "low D4000„" regime (1.8 < £>4000„ < 1.95 for BC03 and 
1.75 < D4000„ < 1.93 for Mastro) and one characteristic of the "high £>4000„" 
regime (1.95 < £>4000„ < 2.05 for BC03 and 1.93 < L>4000„ < 2.05 for Mastro). 
We demonstrate that, at the condition of studying the model separately in the two 
regimes, the linear approximation: 

D4000„(Z) = A(Z) ■ age + B(Z) (2) 

is valid and accurate at fixed metallicity, having in the case of BC03 models correlation 
coefficients always > 0.996 with a mean value of 0.9987 ± 0.0004 for the "low D4000„ " 
regime and always > 0.994 with a mean value of 0.9976±0.0007 for the "high £>4000 n " 
regime; in the case of Mastro models the correlation coefficients are always > 0.997 
with a mean value of 0.9984 ±0.0003 for the "low £>4000„ " regime and always > 0.993 
with a mean value of 0.9986 ± 0.0005 for the "high D4000„" regime. 

The only limitation of this approximation is that it requires a direct estimate of 
the metallicity of the sample. However, as we will show later, this method has been 
proven to be robust against the metallicity evaluation, since the analysis of subsample 
with a median metallicity completely different give results in full agreement. We 
decided to define for each metallicity a mean slope < A >, averaging between the 
slopes with same metallicity and different SFHs. The mean values of the slopes 
are reported in Table [2J The use of an average slope instead of the single values 
is motivated by the fact that the dependence on the SFH is much less significant than 
the dependence on age and metallicity (see upper panels of Fig. [4j, as testified by the 
small errors associated to the mean values (see Table [2]). In the lower panels of Fig. 
|4] we plot in gray the models, as explained before, while the colored lines represent 
the average linear fit to those models; the comparison between the colored and the 
grey lines shows visually the similarity between the average slopes and the slope of 
each single model. The dotted horizontal line represent the separation between the 
"low -D4000„ " and the "high Z?4000 n " regime. From this figure are clearly evident 
the different slopes present in the two range of Z?4000„, for each metallicity. 

Since, as Eq. [2] shows, the slope A represents the conversion parameter between 
the .D4000 n and the age of a galaxy, a fundamental step of our approach is to assign to 
correct A value to the metallicity evaluated in each sample. Therefore we decided to 
interpolate the values reported in Tab. [2] with their errors with a quadratic function, 
obtaining a < A >-metallicity relation; in this way it is possible to find the correct 









low D4000„ 


high D4000„ 


BC03 
BC03 
BC03 


<A> (Z/Zq -- 
<A> (Z/Z G 
<A> (Z/Z G -- 


= 0.4) 

= 1) 
= 2.5) 


0.02893 ± 0.00002 
0.0601 ± 0.0006 
0.1926 ±0.0113 


0.0159 ±0.0003 
0.0359 ±0.0005 
0.1223 ±0.0013 


Mastro 
Mastro 
Mastro 


<A> (Z/Z e -- 
<A> (Z/Zq 

<A> (Z/Zq 


= 0.5) 

-1) 
= 2) 


0.03136 ±0.0004 
0.0582 ±0.002 
0.135 ±0.0083 


0.02427 ±0.00005 
0.0341 ±0.0005 
0.1257 ±0.004 



Table 2. Mean slopes < A > of the D4000„-age relation (see Eq. for the BC03 
and Mastro models with different metallicity and for the two D4000 n regimes. 
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< A > parameter considering the metallicities of each mass subsamplesusing this 
relation, and they are given in Tab. [TJ 

In our approach we take into account also the effect of averaging the slope A between 
different star formation histories, since the quadratic function has been defined with 
a best fit to the < A >-metallicity relation, where the errors in the parameter 

< A > quantify exactly the difference of the slope between the different star formation 
histories considered. The size of the errors by themselves shows that the dependence of 
the slope on the star formation history is small. We find that the percentage difference 
between models with different SFHs, averaged in the metallicity range spanned by our 
data, is 2.8% in the "low L>4000„" regime and 3.1% in the "high £>4000„" regime for 
BC03 models and 12.4% in the "low D4000 n " regime and 3.9% in the "high D4000„" 
regime for Mastro models. We decided to take into account this difference in our 
analysis, adding these errors in the cosmological parameter evaluation. 

This analysis has been carried out also with the new Mastro models to study 
the dependence of the cosmological parameter evaluation on a particular library of 
models. By comparing the plots from Fig. [4] it can be found that the -D4000„-age 
relations have a different normalization, yielding the Mastro models to younger ages 
on average. However, even with different absolute Z?4000„ values, the study of the 
slope of the relations gives values in good agreement with the ones obtained from the 
BC03 models. 

This results is of key importance in our approach. On one hand it is fundamental to 
remind that our analysis is based only on the relative change of D4000„, and not on its 
absolute value. Thus, the variation in the normalization between the different models 
is irrelevant. On the other hand, the fact that the slope of the D4000„-age relation 
remains almost unchanged passing from BC03 to Mastro models demonstrate the 
robustness of our method against the choice of different stellar population synthesis 
models. In Section 15.41 we will quantify what is the effective consequence of this 
difference in the parameter estimation. 



3.2. A linear model 

The cosmological expansion history is described by the Hubble parameter evolution: 

- -This < 3 > 

The analysis of the differential age evolution of ETGs, that trace the differential 
age evolution of the Universe, fully determine H(z), and therefore the cosmological 
parameters. 

As shown in section [3~T1 there exists, in the range of ZM000„ probed by our data and 
at fixed metallicity Z, a linear relation between the D4000 n and the age of a galaxy; 
therefore, if we use Eq. [21 Eq. [3]can be easily rewritten as a function of the differential 
evolution of the D4000„ with a proper conversion parameter, where this conversion 
parameter is given by the slope A of the Z?4000„-age relation: 

H ^ = -i^^m n (4) 

This approach, provided that the galaxies are passive and evolved passively, is 
promising because: 

• it relies on a direct observable of a galaxy, that can be easily measured from 
spectra; 
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• it depends only on the slope of the £>4000 n -age relation, and not on its overall 
normalization; 

• it depends only marginally on the SFH assumptions; 

• it is not influenced by dust extinction (Hamilton 1985, Balogh et al. 1999). 

Equation U can be reformulated in order to make explicit the dependence on the 
various cosmological parameters in the following way: 

#4ooo„(z) = 4^1^ y 1 ' 2 { fi ™ + Qde [y-^o+ Wa ) e3wa{y -i)j y 1 ' 2 dy 

+ constant (5) 
where we have used the standard parametrization w(z) = wo + {jr^)w a . 

4. Deriving the £>4000„-redshift relation 

For each mass subsample, we evaluated the median Z34000 n in narrow redshift bins 
(Az = 0.01); the associated errors are standard errors on the medianJJ Our results are 
shown in Fig. [SJ In the lower mass range, we limit our analysis to z < 0.24, because, 
given the magnitude limit of our sample, above z « 0.24 we have only N ga i = 9. 
It is impressive that for each mass subsample we find a clear Z?4000 n -z relation. Since 
we have proven in Sect. 12.21 that the metallicity evolution with redshift is negligible 
within each mass subsample, this directly witnesses the differential age evolution of 
ETGs, making evident that galaxies at low redshift have always higher break with 
respect to galaxies at high redshift. Moreover Fig. [5] provide also good evidence of 
mass-downsizing, showing that more massive galaxies present at each redshift higher 
Z?4000 n breaks with respect to less massive ones, hence having higher ages, since the 
difference in metallicity is small. 

Comparing our Z?4000„ values shown in Fig. [5] with theoretical values, we find 
that with these metallicities the relative age evolution of the ETGs is w 1.2 Gyrs in 
the range 0.15 < z < 0.3, fully compatible with the theoretical age evolution expected 
in a ACDM Universe within this redshift range. 

The absolute ages of these ETGs go on average from « 5 — 6 Gyrs in the lowest mass 
bin to w 6 — 7 Gyrs in the highest mass bin, so it seems that there are some problems of 
normalization between the ages obtained from the D4000„ and the expected absolute 
ages, that are around ~ 8 — 10 Gyrs for ETGs in this redshift range. However on the 
one hand, it is fundamental to remember again that our approach relies completely on 
the differential Z?4000„ evolution of ETGs, so in the cosmological parameter evaluation 
we are insensitive to the overall normalization. On the other hand, it is also important 
to stress that we are using median values, so we are averaging between different 
galaxies ages; if we consider only the upper envelope of the D4000 n distribution, 
which represents the oldest population of ETGs, we retrieve again absolute ages of 
?» 8 Gyrs. 

Another crucial point in favor of using our approach is that the use of this linear 
conversion from Z)4000„ to the age of a galaxy keeps by definition unchanged the slope 
of the redshift dependence, allowing a better estimate of cosmological parameters. 

f The error on the median are evaluated as the median absolute deviation/sqrt(N), where the "median 
absolute deviation" (MAD) MAD = 1.482 * median (|D4000„ - median(D4:000„)\). (see Hoaglin et 
al. 1983) 
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Figure 5. D4000 n -redshift relation for ETGs in different mass sub-samples. 

As stressed before, to apply our technique it is fundamental to obtain the 
correct A parameter, that is the conversion factor from the D4000„ to the age 
of a galaxy. Therefore we need to associate the proper metallicity to each mass 
subsample. Since, as already underlined, we verified that the metallicity evolution with 
redshift is negligible within each subsample, we studied the medians of the metallicity 
distribution, with their standard errors. As predicted by Gallazzi et al. (2006), we 
find that all ETGs have slightly oversolar metallicity values, increasing as a function 
of mass. 

We use the < A >-metallicity relation defined in Sect. 13.11 to find the proper < A > 
parameter corresponding to the metallicity of each mass subsample, giving the median 
metallicity reported in Table [T] We decided to use the values corresponding to the 
"low L»4000„" regime for the lower mass bins (10.75 < log(M/M Q ) < 11, 11 < 
log(M/M G ) < 11.25, and 11.25 < log(M/M Q ) < 11.5), and the values corresponding 
to the "high D4000„" regime for the highest mass bin (log(M/M Q ) > 11.5). 

5. Constraints on the cosmological parameters 

The Z)4000„-z relations obtained from our selection of ETGs extracted from the 
SDSS-DR4 are then used to set constraints on the Hubble constant Ho and on the 
dark energy equation of state parameter w, assuming w = constant (w = wq, w a = 0). 

The most recent evaluation of the Hubble parameter has been derived by Riess 
et al. (2009) by analyzing the magnitude-redshift relation of 240 low-z Type la 
supernovae at z < 0.1. The absolute magnitudes of supernovae are calibrated 
using new observations from HST of 240 Cepheid variables in six local Type la 
supernovae host galaxies and the maser galaxy NGC 4258. This refurbished distance 
ladder based on extensive use of differential measurements allowed them to obtain 
Ho = 74.2 ± 3.6 km s~ r Mpc , including both statistical and systematic errors. 
This measurement has been proved to be extremely useful also for a new determination 
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Figure 6. D4000 n -redshift relation (upper panels) and dz/dDA000 n relation for 
different theoretical lines, in the upper panels the black points are the relation 
obtained for the 11 < logM / 'Mq < 11.25 mass subsample. The blue lines 
represent the best-fit to the data, using !!m = 0.27, SIde = 0.73, Ho = 74.2 
in the left panels and Qm = 0.27 and Q,de = 0.73 in the right panels. The green 
and red lines show respectively the curves representing the 68% and 90% errors 
to the best-fit. 



of the dark energy equation of state parameter obtained by the seven-year WMAP 
observation (Komatsu et al. 2010): by combining the WMAP data with the 
latest distance measurements from the Baryon Acoustic Oscillations (BAO) in the 
distribution of galaxies (Percival et al. 2010) and the Hubble constant (H ) 
measurement (Riess et al. 2009) it was possible to set the constraint w = —1.10 ± 0.14 
(68% CL). 

The most recent estimate of w comes from the analysis of a sample of 557 SNe called 
Union2 compilation, obtained combining the Union compilation with the HST-light 
curves and spectra of six new high redshift SNe la (Amanullah et al. 2010); their 
analysis gives results in full agreement with WMAP estimate, with a w = — 0.997±0.08 
for a flat ACDM Universe and w = -1.038±0.09 for a ACDM Universe with curvature, 
both assuming w ^constant. 

In the next sections we will show the strength of our approach in determining both 
Ho and w, by analyzing individually the different mass subsample and then performing 
a joint analysis, unifying the constraints derived from the single measurements. 
As an example, we show in Figure |6] in black the D4000„ — z relation for the mass bin 
11 < logM/Mo < 11.25. We performed a best-fit to this relation, in order to obtain 
w (in the left plot) and Hq (in the right plot). The blue curves show the best- fits 
to the data; we also plot in green and red the curves representing the 68% and 90% 
errors to the best-fits, respectively. In the upper-left panel, the best-fit model has been 
evaluated at fixed = 0.27, SIde = 0.73 (as obtained by WMAP 7-years analysis), 
Hq = 74.2 (as obtained by Riess et al. 2009), to set constraint on w = constant; in 
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Figure 7. \ 2 °f the Hq (left panel) and of the w (right panel) evaluation for 
ETGs in different mass subsamples using BC03 models. 



the upper-right panel we set w = — 1, 0m = 0.27, SIde = 0.73, and search for the 
best-fit Ho value. The size of the errors of our data makes clear the strength of the 
differential age approach to set constraints on the cosmological parameters. In the 
lower panels we show the differential evolution of the D4000 n , dz / dDA000 n , showing 
how much these relations are sensitive to the variation of Hq and w. 

5.1. Statistical methods 

We fit the D4000 n — z relations shown in Fig. [5] with the formula given in Eq. [4j 
using a standard % 2 approach. The theoretical £>4000„(z) relation is a function of 
the A parameter and of various cosmological parameters (H Q , Cl m , £Ide, w, w a ). 
Therefore we will have X 2 ( a ii a 2i ■•■ 1 Q ; Af); where are the free parameter of the fit. 
The associated likelihood is given by C(ai, a.?, ckat) oc e~ (x 2 («i. a 2)--,a-w))/2_ Being 
interested only in the evaluation of one parameter, we decided to marginalize this 
likelihood over the parameters we are not interested, and we obtain: 

f X 2 ( a l Q 2 a N) 

C(a.i) = / e 5 dai...don-\donj r i...daN (6) 

In the following sections we will explain the adopted priors and assumptions to 
estimate the Hubble constant Hq and the dark energy equation of state parameter 
w. 

5.2. The Hubble constant 

To estimate H , we fitted the D4000„-z relations given in Fig. [S]with Eq. |U fixing 
as a prior only a flat AC DM cosmology (O m + Ude = 1, w = —1 and w a = 0), 
letting ftoE free to vary in the range provided by the latest WMAP 7- years analysis 
(0.7 < SIde < 0.76) and the slope A free to vary within the values associated to the 
metallicity range found in the data within la error (see Table [T]), considering also the 
spread due to the different SFHs considered (see the second point of Sect. 15. 4j) . 
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As explained in Sect. 15.11 the best-fit has been evaluated with a standard x 2 
analysis, marginalizing the obtained likelihood over those parameters which we are 
not interested to constrain, that are SIde and A. We end up with four likelihoods of 
the Hubble constant for the mass bins we created. The estimates obtained from the 
higher populated samples (11 < log(M/M Q ) < 11.25 and 11.25 < log(M/M Q ) < 11.5) 
show a perfect agreement, with comparable errors; the analysis of the other two mass 
subsamples give slightly different results, but still within the la errors. However, it is 
also important to remember that the number of galaxies of these subsamples is quite 
lower with respect to the bin 11 < log{M / M Q ) < 11.25 (respectively 36% and 72% 
less), and therefore that the associated errors are larger (see Tab. [3]). 

Since those likelihoods are independent from each other, to reduce the statistical 
error and to obtain a better estimate on H we decided to join them, obtaining 
H = 72.6 ± 2.9 kmMpc _1 s _1 . In Fig. [7] are shown the x 2 for each mass bin and the 
joint x 2 , and in Tableware reported the corresponding values of Hq and their errors. 



5. 3. The dark energy equation of state 

To estimate the value of w, assuming w — constant, we fitted the Z?4000 n -z relations 
using as priors the VI de value as found by the latest WMAP 7-years analysis 
considering a wCDM flat Universe (Larson et al. 2010, f2 m + Ude — 1, Wa = 
and 0.645 < VIde < 0.835) and the estimate of Hq done by Riess et al. (2009) 
(70.6 < Hq < 77.8), the slope A free to vary within the values associated to the 
metallicity range found in the data within la error (see Table Q]), considering also the 
spread due to the different SFHs considered (see the second point of Sect. I5.4j) . 
The best-fit has been evaluated with a standard \ 2 analysis, marginalizing afterwards 
over H , VLde and A. We end up with four likelihoods of w for the mass bins we 
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Table 3. Hubble constant (in unit of kmMpc J s 1 ) and relative errors using 
BC03 models. 
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Table 4. Dark energy equation of state parameter w and relative errors using 
BC03 models. 
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created. All the estimates obtained on w are in agreement within la error; only the 
value obtained in the highest mass bin shows a slight discrepancy, even if compatible 
with the other estimates within la. However, a similar discussion as the one introduced 
in the Hq estimate has to be done for this discrepancy, considering that it is the bin 
with the highest errors and the lowest galaxy number. As done in the Hq analysis, 
since those likelihood are independent one from the other we performed a joint analysis 
by multiplying them, to reduce the statistical error. In Fig. |9] are shown the x 2 for 
each mass bin and the joint x 2 , and in Table |4] are reported the corresponding values 
of w and their errors. From the joint analysis we obtained w = —1 ± 0.2. 

5.4- Study of the systematics 

In our approach, we studied the effect of different possible systematics from which our 
analysis may suffer. In particular, the most critical ones are the dependence of our 
results on the estimate of the absolute age of a galaxy from the D4000„, on the star 
formation histories considered, on the stellar population synthesis models assumed, on 
the metallicity evaluation and on the not perfect assumption of the passive evolution 
of our sample. In the following we address separately all these sources of uncertainty 

• Dependence on the absolute age evaluation. Since our approach relies totally on 
the study of the differential age (or Z?4000„) evolution to constrain cosmological 
parameters, our approach is insensitive to a possible bias in the evaluation of the 
absolute age of a galaxy. 

On the contrary, we strongly rely on the fact that our technique does not smear in 
some way the slope of the age-redshift evolution, i.e. that the differential evolution 
is preserved unmodified. This approach, by definition, keeps the differential 
evolution unchanged, since we use a simple linear conversion between the Z?4000„ 
and the age of the galaxy, justifying in Sect. 13. li the robustness of this assumption. 




Figure 8. x 2 °f the Ho (left panel) and of the w (right panel) evaluation for 
ETGs in different mass subsamples using Mastro models. 
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Table 6. Dark energy equation of state parameter w and relative errors using 
Mastro models. 



• Dependence on the star formation histories considered. In Sect. 13. II we justify the 
choice of using delayed exponentially declining star formation histories with low 
r's by analyzing the SED-fit tau distributions of our sample (see Fig. [5]). However, 
the study of models with different r's introduces a spread in the conversion factor 
A between the ZM000„ and the age of a galaxy, since we averaged the values found 
with the different SFHs at fixed metallicity. The percentage spread in the slope, 
averaged in our metallicity range, is 3.1% and 2.8% for BC03 models respectively 
in the "high D4000" regime and 13.4% in the "low D4000" regime), and 3.9% 
and 12.4% for Mastro models, and the effect is even smaller (almost negligible) 
in the estimated cosmological parameters. Even if this spread is small, we take 
in consideration this additional error in our analysis. 

• Dependence on the choice of the stellar population synthesis model. In Sect. I3.1l wc 
perform the analysis of the slopes of the D4000„-age relation both using BC03 and 
Mastro models. As pointed out before, these two models are completely different, 
being the Mastro models based on the latest MILES models and the BC03 models 
on the STELIB ones. In Tab. [2]we show the values of the slopes in the two D4000„ 
regimes for the two models. To quantify the difference between these two model, 
since we are interested in our particular metallicity range shown in Tab. [T] we 
interpolate the values of the slope < A >-mctallicity relation in Tab. [5] for both 
models with a quadratic function. We find that the mean percentage difference 
between the slope obtained from the two models in our metallicity range is 2.8% 
in the "low D4000" regime and 13.4% in the "high D4000" regime. 

We also evaluated the effect of this percentage difference between the models by 
performing the same analysis described in Sect. 15.21 with Mastro models. The 
values of Hq and w obtained in the different mass bins are shown in Tab. [5] and 
Tab. El and the corresponding x 2 ar< 3 shown in Fig. [8j After marginalizing and 
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H 


68% errors 


90% errors 


# galaxies 




km Mpc _1 s _1 


km Mpc _1 s _1 


km Mpc _1 s _1 




11 < log{M/M & ) < 11.25 


72.2 


+4J/-4.2 


+8.3/-6J 


5429 


"low Z" 


71.7 


+7.7/-6.1 


+ 13.5/-9.6 


2717 


"high Z" 


70.8 


+5.3/-4.S 


+9.1/-7.2 


2712 


EWo(Ha) < 0.5A 


70.4 


+5.8/-4.8 


+10/-7.7 


3245 



Table 7. Hubble constant (in unit of kmMpc -1 s~ 1 ) and relative errors for 
11 < log(M / Mq) < 11.25 sample, for the same mass bin sample divided below the 
median mctallicity, above the median metallicity and selected with EWo(Hct) < 
0.5A, using BC03 models. 





w 


68% errors 


90% errors 


# galaxies 


11 < log(M/M Q ) < 11.25 


-1 


+0.3/-0.5 


+0.6/-0.8 


5429 


"low Z" 


-1.1 


+0.5/-0.6 


+0.9/-1.1 


2717 


"high Z" 


-1.2 


+0.4/-0.5 


+0.7/-0.9 


2712 


EW {Ha) < 0.5 A 


-1.2 


+0.4/-0.5 


+0.7/-1 


3245 



Table 8. Dark energy equation of state w and relative errors for 
11 < log(M I Mq) < 11.25 sample, for the same mass bin sample divided below the 
median metallicity, above the median metallicity and selected with EWo(Ha) < 
0.5A, using BC03 models. 



doing the joint analysis combining the informations of all the mass bins, we find 
using Mastro models H = 73.4 ± 3.5 kmMpc"^" 1 and w = -0.9 ± 0.3. We find 
therefore no significant difference between BC03 and Mastro models, obtaining 
values for both Hq and w compatible within the la errors. 

• Dependence on metallicity and on star formation. An important test of this 
analysis is to verify the robustness of our results against the choice of the 
value of the mctallicity, since it is the primary parameter to obtain the correct 
slope of the Z)4000, i -age relation. To check that our analysis is not biased 
in that sense, we studied the mass bin with the highest number of galaxies 
(11 < log(M /Mq) < 11.25), dividing it into two sub-sample sampling different 
metallicity ranges; we decide to divide the sample below its median metallicity 
{"low metallicity" sample) and above its median metallicity {"high metallicity" 
sample), to keep the number size of the resulting subsamples similar. In this way 
we end up with two subsample with significantly different values of metallicity 
{Z/Zq = 0.92 and Z/Z Q = 1.38, respectively). 

In Fig. [§] (upper panel) it is shown the Z)4000, i -z relation of the parent 
sample (in black) and of the "low metallicity" and "high metallicity" subsamples 
(respectively in blue and red). We evaluated the median metallicity of the two 
subsamples, to obtain for each one the correct A parameter. The corresponding 
Z?4000„-redshift relations result rather different between the "low metallicity" 
and the "high metallicity" sample; however given the different metallicity and 
therefore conversion parameter A, the fit to Hq and w gives results in perfect 
agreement, as it is possible to see from the x 2 in the lower panels of Fig. [S] and 
from the values reported in Table [7] and Table [5J 

Furthermore, from the analysis of the emission lines of our ETGs, we found 
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Figure 9. In the left panel is shown the D4000-rcdshift relation for ETGs with 
11 < log(M I Mq) < 11.25 (in black), for the same mass bin sample divided below 
the median metallicity (in blue) and above the median metallicity (in red), and 
for the same mass bin selected with EWo(Ha) < 0.5 A (in green). In the right 
panel there are the corresponding x 2 ■ 



that, in spite of the rather strict selection criterium applied, we still have a tail 
in the EW (Ha) and EW {[0 1 7]A37 '27). On the one hand, as pointed out in 
Sect. 12.21 these are quite low levels of equivalent widths, and the distributions 
of the eClass parameter and the specific star formation rate obtained from the 
SED fitting for the galaxies with a significant detection are characteristic of a 
passive population. On the other hand, this may be an indication that our 
sample contains a certain number of galaxies with undergoing star formation 
or AGN activity. Our approach relies strongly on the assumption that we are 
selecting a sample of passively evolving galaxies, that traces uniformly the age 
evolution of the Universe as a function of redshift. To test if there is a such kind 
of bias in our results, we decided to apply to the same mass subsample chosen 
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Source of uncertainty 


Effect 


Impact on 


Impact on 


uncertainty 




Ha 


w 


Absolute age evaluation 


possible underestimate of the 
absolute ages using D4000„ 


none 


none 


SFHs assumption 


spread in the conversion factor A 


< 0.01% 


< 0.1% 


Stellar population 


different D4000„-age relation 


1.1% 


10% 


synthesis model used 


within models 






Metallicity estimate 


if incorrect, may bias the A parameter 


1.4% 


15% 


non-perfect approximation 


possible bias in the D4000„-z evolution 


2.5% 


20% 


of passive evolution 








Total systematic error 




3.1% 


27% 



Table 9. Summary of the sources of uncertainty and their impact (in percentage) 
on the cosmological parameters. 



before (11 < log(M /M©) < 11.25) an even stricter selection, considering only 
those galaxies for which EWo(Ha) < 0.5 A (that is the median of the EWo(Ha) 
distribution). We evaluated the metallicity of this sample, and redid the analysis 
of the Hubble constant and dark energy equation of state parameter, and the 
results are shown in Fig. |9] (green lines), and Tab. [7] and Tab. HI We find that 
the value obtained from this analysis well agrees with the value obtained not 
cutting the sample in EWq (Ha) . So we conclude that our results are robust also 
against the presence of galaxies with small values of EWo(Ha). 

6. Total systematic uncertainty and results 

In Sect. 15.41 we quantified the impact of systematics on the estimate of both the 
Hubble parameter Hq and the dark energy equation of state parameter w (assuming 
w = constant). All these source of uncertainties, along with their effect, are 
summarized in Tab. [9j to evaluate the effect of systematics, most of the tests 
have been performed using the mass sample with the highest number of galaxies 
(11 < log(M /Mq) < 11.25). The total estimated effect has been evaluated by 
summing in quadrature the errors, and it is also reported in Tab. |H1 

In conclusion, we find that our estimated value for Hq, including both 
statistic and systematic errors, is Hq = 72.6 ± 2.9(stat) ± 2.3(syst) kmMpc _1 s _1 ; 
our estimated value for w, including both statistic and systematic errors, is 
w = -1 ± 0.2(stat) ± 0.3{syst). 

7. Summary and conclusions 

We have developed a new technique to obtain the expansion rate of the Universe using 
the "cosmic chronometer" technique with the aim of minimizing the dependence on 
systematics. To this extent, we have shown that the D4000„ feature at fixed metallicity 
correlates linearly with age, for the range of ages of interest. We have studied this 
feature using theoretical synthetic stellar population models and shown that it is 
robust to choice of metallicity, star formation history and different stellar population 
models. We also evaluated, within this linear approximation, the theoretical model 
Z34000„(z) as a function of the various cosmological parameters. 
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We obtained a sample of ETGs from the SDSS survey, using both photometric 
and spectroscopic information to select the most passive and massive galaxies 
(log(M/M Q ) > 10.75). We divide our sample in small mass bins (0.25 dex), in order 
to avoid possible bias from mass downsizing and to have homogeneous redshift of 
formation and metallicity across the entire redshift range. The Z)4000 n — z relations 
show, in each subsample, a clear redshift evolution, with galaxies at low redshift having 
always a higher D4000„, and hence higher ages. Moreover the four subsamples give an 
evident proof of mass-downsizing, since at each redshift galaxies with higher masses 
have always an higher D4000„ than galaxies at lower masses. 

We evaluate the Hubble constant H and the dark energy equation of state param- 
eter w by fitting the D4000„ — z relations with the theoretical model defined before. 
We find from the joint analysis that H = 72.6 ± 2.9(stat) ± 2.3(syst) kniMpc^s -1 
and that w = — 1 ± Q.2(stat) ± 0.3(syst), assuming a constant w. The value we ob- 
tained for w is in agreement with the value found from the WMAP 7-years analysis 
and with the one obtained from Amanullah et al. (2010); the value of H n is within la- 
error with the one obtained by Riess et al. (2009), and we reach a comparable level 
of precision. We have also shown that the universe is accelerating, even if our sample 
only reaches up to z ~ 0.3. 

We demonstrate that our results are robust to the choice of galaxy mass, since 
the results in the single mass subsample are all within the ler errors. Moreover, 
we study the effect of the dependence of our results on different systematics, i.e. 
the estimate of the absolute age of a galaxy from the D4000„, the star formation 
histories considered, the stellar population synthesis models assumed, the metallicity 
evaluation and the not perfect assumption of the passive evolution of our sample. 
We find a small dependence on the SFHs studied, that we took into account in the 
cosmological parameters evaluation. By studying our sample with different synthetic 
stellar populations models (BC03 and Mastro), in different metallicity regimes and 
for different selection of EWo(Ha), we show that our analysis is also robust to the 
metallicity evaluation and to the non-perfect approximation of passive evolution of 
our sample, and that there is no significant difference when using different models for 
studying the D4000„-age relation. 

In forthcoming papers we will apply our method to higher redshift samples, to 
determine the expansion history of the Universe, i.e. H(z), up to z ~ 1. This method 
is promising also in view of the future massive spectroscopic surveys expected in the 
next decades from the ground (e.g. SDSS-III BOSS, BigBOSS) and from space (e.g. 
ESA Euclid, NASA WFIRST). 
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